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Abstract 

We have used data from the TeV 7-ray flare associated with the active 
galaxy Markarian 421 observed on 15 May 1996 to place bounds on the 
possible energy-dependence of the speed of light in the context of an 
effective quantum gravitational energy scale. The possibility of an ob- 
servable time dispersion in high energy radiation has recently received 
attention in the literature, with some suggestions that the relevant en- 
ergy scale could be less than the Planck mass and perhaps as low as 
10 16 GeV. The limits derived here indicate this energy scale to be in 
excess of 4 x 10 16 GeV at the 95% confidence level. To the best of our 
knowledge, this constitutes the first convincing limit on such phenom- 
ena in this energy regime. 
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It has recently been pointed out that many quantum gravity scenarios may 
result in an observable time dispersion for high energy radiation originating at 
large distances from the Earth QTJ This would result from an effective 

energy-dependence to the velocity of light in vaccum owing to propagation 
through a gravitational medium containing quantum fluctuations on distance 
scales near the Planck length, Lp ~ 10 _33 cm, with timescales on the order of 
1/Ep, where Ep is the Planck Mass (~ 10 19 GeV). In particular, it has been 
indicated |l[ that different approaches to quantum gravity lead to a similar 
description of the first-order effects of such a time dispersion: 

A***/-- (1) 

Eqg c 

where At is the time delay relative to the standard energy-independent speed 
of light, c; £ is a model-dependent factor of order 1; E is the energy of the 
observed radiation; Eq G is the assumed energy scale for quantum gravitational 
effects which can couple to electromagnetic radiation; and L is the distance 
over which the radiation has propagated. While Eqg is generally assumed to 
be on the order of Ep, recent work within the context of string theory suggests 
that the onset of noticeable quantum gravitational effects may correspond to 
a characteristic energy scale smaller than the Planck mass and perhaps as low 
as 10 16 GeV Q. Thus, any experimental probe of such scales or higher would 
be of great interest. 

In a recent paper |3J it was suggested that 7-ray bursts (GRBs) could 
provide a natural way to test such predictions owing to the short duration, high 
energies and the apparent cosmological origin of at least some of these bursts. 
Based on current data, these authors indicate that if 1) time structure on the 
scale of 0.01 seconds or smaller can be established for energies ~ 200keV and 2) 
an association of such a burst can be made with an object possessing a redshift 
of order 1, energy scales of Eqg ~ 10 16 GeV could be probed. Unfortunately, 
establishing the distance of any particular GRB from earth has proven to be 
non-trivial, with only a handful positively associated with optical counterparts. 
Also, some of the highest energies seen from GRBs are associated with an 
"afterglow" which seems to occur over much longer timescales than the initial 
burst. However, more stringent and robust limits to Eqg can already be set 
based instead on the rapidly rising TeV flares seen to occur in active galaxies. 

The Whipple Observatory 7-ray telescope, located in Arizona, detects the 
Cerenkov light generated by electromagnetic cascades resulting from the inter- 
action of high-energy 7-rays in the atmosphere. Images taken of such cascades 
are used to discriminate backgrounds and derive energies of the primary 7- 
rays in the regime above ~250 GeV. To date, three extragalactic sources, all 
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active galaxies of the blazar class, have been identified as emitters of TeV ra- 
diation Two of these, Markarian 421 and Markarian 501, produce 
particularly strong emission with energy spectra approximated by an ~ E~ 2 5 
power law (although Markarian 501 shows evidence for additional curvature) 
between energies of 300 GeV and 10 TeV |§. These same sources have also 
exhibited dramatic changes in flux level on timescales ranging from minutes to 
days. On several occasions, such variations have been simultaneously studied 
and correlated with x-ray, UV and optical measurements |§ |L0| . 

The most rapid flare observed thus far was seen from Markarian 421 on 
15 May 1996 [O]. This data is shown in figure 1, where the excess rate 
of 7-ray selected events above a threshold of 350 GeV is binned in intervals 
of 280 seconds duration, as it appeared in the original publication of this 
observation. To avoid confusion (and potential bias), we will retain this same 
binning for the current analysis. The doubling time of the flare is less than 15 
minutes, although variability is apparent on the scale of the binning at the 99% 
confidence level. Because of the rapidly falling energy spectrum, the 7-ray data 
is dominated by events near the triggering threshold. Thus, the peak of the 
flare is almost entirely defined by events with 7-ray energies less than 1 TeV, 
as shown at the top of figure 2. In the 280 second interval corresponding to 
this peak, 4 events with 7-ray energies in excess of 2 TeV have been identified, 
whereas no such events could be identified in either of the 2 adjacent intervals 
(figure 2, bottom). The probability for this to occur by chance is less than 
5%. Hence, at the 95% confidence level, emission above 2 TeV appears to 
keep in step with emission below 1 TeV for variability timescales less than 280 
seconds. The redshift of Markarian 421 is 0.031, which translates to 1.1 x 10 16 
light-seconds for an assumed Hubble constant of 85 km/s/Mpc. From equation 
1, this leads to lower bound on Eqg/£ of 4 x 10 16 GeV. 

We note that an earlier limit on the energy- dependence of the speed of 
light, which would be more restrictive than that given here, had been derived 
from the possible ultra-high-energy detection of anomalous pulsed emission 
from Her X-l in 1986 |Tj|. However, more recent analyses and the lack of 
further such detections suggests that the interpretation of that observation as 



a statistical fluctuation is not an unreasonable one ||13|| . We therefore believe 
that the limit presented in this paper represents the most credible and stringent 
bound thus far obtained. 

The next generation of proposed ground-based instruments, such as VER- 
ITAS and HESS, will feature multi-telescope systems with much improved 
sensitivity, energy coverage and resolution, along with the ability to track can- 
didate sources of flares more continuously using dedicated telescopes. This 
will allow for both a more detailed study of the time structure of currently 
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known TeV sources and the prospect of discovering and studying more distant 
objects. It is therefore reasonable to expect to probe Eqq to even higher ener- 
gies in the near future from further studies of TeV flares. As has already been 
pointed out , the distinctive dependence of the shortest observable variabil- 
ity timescale on both energy and source distance for quantum gravitational 
dispersion should allow source-specific effects to be distinguished. Thus, fu- 
ture TeV studies could conceivably provide convincing evidence for quantum 
gravity, particularly if the resulting time-dispersion effects are associated with 
characteristic energy scales less than the Planck mass. We hope that this 
prospect, in addition to the bounds derived here, will encourage more detailed 
predictions of such phenomena to be calculated in the context of specific quan- 
tum gravity frameworks. 

We would like to thank Subir Sarkar at Oxford for several useful discussions 
on this topic. This work has been supported in part by PPARC, Forbairt, the 
US Department of Energy, NASA and the Smithsonian Institution. 
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Figure 1: TeV 7-ray flare from Markarian 421 observed on 15 May 1996 by the 
Whipple 7-ray observatory. The rate of excess 7-ray selected events is binned 
in intervals of 280 seconds, (taken from reference [H|) 
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Figure 2: Total number of 7-ray selected events occurring in each 280 second 
interval near the peak of the 15 May 1996 flare from Markarian 421. The top 
plot consists of events with 7-ray energies less than 1 TeV, whereas the bottom 



plot is for energies greater than 2 TeV 
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